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Abstract 

We study reflection and transmission of waves in a random tight-binding system with absorption 
or gain for weak disorder, using a scattering matrix formalism. Our aim is to discuss analytically the 
effects of absorption or gain on the statistics of wave transport. Treating the effects of absorption 
or gain exactly in the limit of no disorder, allows us to identify short- and long lengths regimes 
relative to absorption- or gain lengths, where the effects of absorption/gain on statistical properties 
are essentially different. In the long-lengths regime we find that a weak absorption or a weak 
gain induce identical statistical corrections in the inverse localization length, but lead to different 
corrections in the mean reflection coefficient. In contrast, a strong absorption or a strong gain 
strongly suppress the effect of disorder in identical ways (to leading order), both in the localization 
length and in the mean reflection coefficient. 

PACS numbers: 72.15.Rn,42.25.Dd, 42.50.Nn 



I. INTRODUCTION 



In this paper we study analytically the coherent reflection and transmission of waves in an 
active one-dimensional disordered system which either absorbs or amplifies the waves. Our 
model is the familiar single-band tight-binding model with random site energies (Anderson 
model) including additional fixed positive or negative imaginary parts describing absorption 
or amplification. As is well-known, the introduction of the imaginary potential destroys the 
time-reversal symmetry of the system. 

The electronic model with absorption may describe annihilation of electrons via electron- 
hole recombinations acting as a complex optical potential in a nearly compensated semi- 
conductor. The amplification model is, of course, meaningless for electrons whose fermionic 
character forbids the presence of more than one electron at a given spatial location. 

On the other hand, the tight-binding model with absorption or amplification due to 
stimulated emission may be appropriate for describing the localization of light waves in 
active photonic band-gap crystals, characterized by a periodic variation of the dielectric 
constant. In particular, the interplay of the (phase coherent) amplification of light waves 
with the process of coherent scattering by random inhomogeneities leading to localization^ 
is of current interest for random lasers 3 -. 

A considerable amount of theoretical work related to the statistics of the transmit- 
tance and of the reflectance of random systems with absorption or gain has already been 
published ^ 5161 'i 8 ? 9 ! 10 ! 11 ! 12 ! 13 ! 14 ! 15 . We feel, however, that important aspects of the effects of 
absorption or gain on the transmission properties of the random system studied below, have 
not received sufficient attention in previous work. 

We refer, in particular, to the form of the localization length £ for large lengths L of a 
random chain. For weak disorder this is found to be given by^ 

7 = -r + T> u = aor 9 , (i) 

? I'll ^0 

where £o is the localization length of the system in the absence of absorption or gain (am- 
plification), l a and l g are the absorption and gain lengths of a perfect system, respectively. 
In the tight-binding model, (JTJ) is expected to be obtained by studying the transmittance of 
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the disordered sample described by the Schrodinger equation 



[E - (e„ + ir})]ip n = y? n+ i + (p n -i, n = 1, 2, . . . 



(2) 



where (p m are the wavefunction amplitudes at sites m = 1, 2, . . . , N, of spacing a, of the 
disordered sample of length L = Na. The e m are the random site energies to which one 
adds a fixed non-hermitian term irj describing absorption for rj > and amplification for 
i] < 0. The energies E, e m and 77 are in units of a fixed nearest-neighbour hopping energy 
V . The random chain is connected at both ends to semi-infinite perfect leads (e m = 0) 
with rj = 0, whose sites are positioned at m = 0, —1, —2, . . . and m = N + 1, N + 2, . . ., 
respectively. For the perfect tight-binding chain with absorption or gain, Datta^ has derived 



while numerical studies of the random tight-binding model by Gupta, Joshi and Jayannavar— 
and by Jiang and Soukoulis 13 support Eq. (JIJ for small 77 of either sign, with l/£o given 
by the familiar Thouless expression for weak site-energy disorder. A feature of (|1J) which is 
generally regarded as paradoxical is the fact that it leads to suppression of transmittance 
for large L, as in the case of absorption. What is more, according to equation (JIJ, when 
disorder is present the suppression for amplification occurs at exactly the same rate as 
for absorption 9 . This surprising feature has been an important incentive for developing a 
more comprehensive analytic treatment of the statistics of wave transport in the presence of 
absorption or amplification. Indeed, an important drawback of Eq. ([TJ is that it completely 
ignores effects of absorption or gain on the statitstics of the transmission coefficient from 
which I is obtained. Our aim is to remedy to this defect in the framework of a detailed 
analysis for weak disorder of the tight-binding system defined by Eq. (J2J). 

An important feature of our approach below is an exact treatment of the effect of absorp- 
tion or gain, as done previously by Datta^ for a non-disordered (pure) system. This allows 
us to clearly identify and discuss short- and long lengths regimes relative to the absorption 
(amplification) length in (|TJ). 

The study of transmission and reflection in random one-dimensional media with ab- 
sorption or gain was initiated, and later pursued actively, using invariant imbedding 
equations ^ 5 i 6 i 7 i 8 i 9 i 10 . These equations are coupled non-linear differential equations for the 
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(3) 



a 
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reflection- and transmission amplitudes of plane waves incident at the right of a continuous 
medium occupying the domain < x < L of the x-direction— . As was shown recently^, 
the invariant imbedding equations follow for weak disorder from the long-wavelength con- 
tinuum limit of (J2J) for a disordered chain embedded in an infinite perfect chain. We recall 
that the invariant imbedding equations were originally derived as an exact consequence of 
the Helmholtz equation for the propagation of the electric field in a dielectric medium ^ 18 ! 19 . 
For later discussion, we also recall the important early result for the so-called short-lengths 
localization length^, 

1 1 

7 = V + T , (4) 
? So 

which indicates that at short lengths transmittance increases with L for amplification, if 

The paper is organized as follows. In Sect. II we derive exact expressions for the scattering 
matrix elements (transmission and reflection amplitudes) in terms of transfer matrices for the 
tight-binding equation (JIJ) for weak disorder. In Sect. Ill we discuss our explicit analytic 
results for the averaged logarithmic transmission coefficients (localization lengths) in the 
short- and long lengths regimes. We also discuss the mean logarithmic reflection coefficient 
for large lengths. We recall that the distribution of the reflection coefficient is important in 
the context of random lasers^- Some concluding remarks follow in Sect. IV. 



II. TRANSFER MATRIX ANALYSIS 

We start our analysis by rewriting (J2J) in terms of a transfer matrix for a site n, 

E — e n — ii] 









|, P„=( 











(5) 

The analogous equation for sites in the perfect leads involves the transfer matrix P° obtained 
by letting e n = rj = in(jSJ). We wish to study the scattering (reflection and transmission) 
of (Bloch)plane wave states of the leads by the disordered segment of length L = N (with 
a = 1). For this purpose it is necessary to perform a similarity transformation of (0) to the 
basis of the Bloch wave solutions tp n ~ e ± lkn f or the leads. The eigenvectors of P° are of the 
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form u± = ( I with eigenvalues e ±lk obeying the equation 

E = 2cosk , (6) 

which defines the tight-binding energy band. As usual we choose k positive, < k < it, so 
that e.g. e lkn corresponds to a wave propagating from left to right on the lattice of (0). The 
similarity transformation of P n to the Bloch wave basis is defined by the matrix U = (u + , u_) 
and leads to 

Q„ = t/- 1 p„f? = g° + g^ , (7) 



where 



and 






' e lk e lk \ 




\-e ik -e- ik ) 



Ql=V ,\,Ql = iK \ , (8) 



" , (9) 



2 sin k 2 sin k 

The transfer matrix of the disordered segment of length iV is the product of transfer 
matrices associated with the individual sites: 

N 

Q=l[Qn ■ (io) 

n=l 

We recall now the precise relationship between the transfer matrix elements (Q J = Qij 

\ / ij 

and the reflection and transmission amplitudes r + ~ and t , and r~ + and t ++ for waves 
incident at the left and at tie right of the disordered system, respectively. The reflection 
and transmission amplitudes define the scattering matrix S, which expresses outgoing wave 
amplitudes at te left (O) and at the right (O') of the disordered segment in terms of ingoing 
ones, (/) and 

7 o\ (r + t—\ ( I 
O' I \ t ++ r + - I \ V 



(11) 



The transfer matrix Q, on the other hand, gives the Bloch wave amplitudes at the right end 
of the disordered section in terms of the amplitudes at the left end: 

(12) 





whose transformation to a form analogous to (jllj) yields: 

o i 1 )(' I ■ < 13 > 

\det Q Q12/ ^ 

which leads to the desired expressions of transmission- and reflection amplitudes in terms 




of the transfer matrix elements Qij = 



t~ = — , t ++ = (detQ)r- , (14) 



Q22' 

r + = — , r + = -— . (15) 

V22 V22 

From (fTU|) it follows that the determinant of Q is the product of the determinants of the 
exact transfer matrices Q n associated with the individual sites n of the disordered segment 
of N sites, of lengths iV a. Now, from (0IHJ) we find that det Q n = 1, n — 1, 2, . . . N so that 

det§=l . (16) 

From ()14j) it then follows that 

t ++ =t-^t = ^~ , (17) 

V22 

for any realization of the disorder and for any strength of the imaginary potential. 

As usual, we assume that the random site energies are identically distributed independent 
Gaussian variables with mean zero and correlation 

For weak disorder, we shall expand the matrix functional Q to linear order in the random 
site energies. We note that since the energies of neighboring sites are uncorrelated (Eq. 
(|18)l) the second order terms in the expansion of (jlO)) may be omitted since they will not 
contribute in averages over the disorder at the order Eq. Thus restricting to first order in 
the expansion of we obtain^ 

N 

Q=(^) iV + E(^) m " 1 ^(^) 7V " m = ^ +^ 1 • ( 19 ) 

m=l 
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The transfer matrix product = ( Q Q n J for the medium in the absence of disorder may 
be readily evaluated in closed form. This will allow us to discuss analytically the reflection 
and transmission properties of a perfectly amplifying or absorbing system, as done earlier 
by Datta^ using a slightly different procedure. We write 



(Qn) =V{V- 1 QIV) V~\ V = (v + ,v.) , (20) 
where V is the diagonalizing matrix formed by the eigenvectors 



/ be~ ik \ 




( be~ ik 


( (l-6)e ifc -e i s 




(l~b)e ik -e- 


V 1 ; 




\ 1 



\ 1 / V 1 / 

(with det V = 2l ^ * sin q) of Ql corresponding to eigenvalues e iq and e~ iq } respectively. 
These eigenvalues are defined in terms complex functions cos q and sin q by 



e ±iq = cosg ± i sing, 



cos q = cos k — ib sin k, sin q = y 1 — cos 2 q . (22) 
The explicit expression of ( ) obtained from (j2()H22)l is 



A R 



Ql) =1""" ~"'| , (23) 

C m Dm , 



where (with e = d±) 



A m = — — [(1 — b)d + singm — sing(m — 1)] , (24) 



D m = [(1 — b)d + sing m — sing(m + 1)1 , (25) 

sin g 



„ , , smgm . . 

B m = -bd„—^— , 26 
sin g 



smgm 

C m = b d + ^ — . 27 
smg 
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Note that, as expected, the transfer matrix of the non- disordered system with absorption or 
gain defined by (j24H27|) obeys 



A N D N - B N C N = cos 2 qN + sin 2 qN = 1 , (28) 
and from (|15j) and f)26M27j) it follows that in this case the reflection coefficients are equal 



\B I 2 

|r 2 | = \r + -\ 2 = \r~ + \ 2 = L^L, e m = 0, m = l,2,...N . (29) 

Finally, by inserting (jHJ) and ()23|) in (fT9")) . we obtain the explicit form of the transfer 
matrix of the disordered section to first order: 




Q\i Q22, 



where 



A? 



m—l 

- B m _ 1 )(d + A N _ m + d-C N ^ m ) , (31) 



m=l 

N 



Q22 — i b m (C m -i — D rn _i)(d + B N _ m + d_D N _ m ) , (32) 



m=l 

N 

Q12 = i b m {A m -i — B m ^i)(d + B N - m + d-D^-m) , (33) 

m=l 
N 

)(d + A N _ m + d„C N „ m ) , (34) 

m=l 

The lowest order effect of the disorder in the mean transmission and reflection coefficients is 
obtained by expanding (|15j) and (|17|) to second order in the disorder, using (J3U|) and noting 
that the first order averages vanish. However in the case of the reflection coefficients, which 
are asymptotically independent of length in the absence of disorder—^, it is apt to focus 
on the simpler averages (ln|r ±T | 2 ). On the other hand, the study of the mean logarithm 
of the transmission coefficient is of special interest since it is related asymptotically to the 
localization length 

1 ,. (In |t ±=t | 2 > 

— = -lim^ — ' ' 1 , 35 

f± 2N y 1 
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which is a self-averaging quantity in the absence of absorption or gain 21 . From ()15jh (|17|) 
and (J3T?|) we obtain successively for the quantities of interest, to second order in the disorder, 




We close this section by demonstrating the equivalence of our results for the transmission- 
and reflection coefficients for the perfect absorbing or amplifying system (e m = 0) and the 
corresponding results obtained earlier by Datta 12 at the band center. For this purpose we 
identify the parameters e iq and e~ iq defined above respectively with the quantities i and 
— i involving the parameter £ introduced by Datta via the substitution sinh£ = |. By 
transforming Eq. (5) o^ for |t| 2 for even N in terms of the variable q we get 



(— zsin Nq + sing cos iVg) 2 ' 
which coincides with the expression obtained by substituting (|2"5j) for k = n/2 and m = N 
in the definition (|17|) of \t |. Similarly, the transformation of Eq. (6) oi^ for odd iV 
yields again (|39*|l obtained from (JT7j) and above. The equations (7) and (8) of Datta^ 
for the reflection coefficient |r| 2 for even and odd N, respectively, reduce similarly to the 
corresponding expression obtained from (|25H26j) and Clearly, the advantage of the 

present treatment is that it condenses distinct expressions for even and odd N in Datta's 
analysis into a single one for any one of the amplitudes coefficients in (I15JI . (|17|). This is 
clearly useful, particularly for handling the more cumbersome general expressions for the 
effect of a weak disorder. 

III. DETAILED RESULTS FOR E = 

For simplicity and as in most previous work for the tight-binding mode l 11 ' 12 ' 13 , we restrict 
the analytical calculations and results in this section to the band center, E = 0(k = vr/2). 
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At the band center the pure system transfer matrix elements (|24H27|) take the simple forms 



Ar, 



iqm 



-iqm 



(40) 



-w_e 



iqm 



+ u+e 



-iqm 



(41) 



B m , = C m = v(e iqm -e~ l ' im ), m=l,2,..., N , (42) 

where 

u± = \( , 1 ± 1 ) , v = - b , (43) 

e ±iq = (±i)(VTT¥^b) . (44) 

For the pure tight-binding system with absorption or gain [e m — 0, m — 1, 2 . . . iV) the 
transmittances and reflectances, for both directions of incidence, are given exactly for any 
band energy E and for any length L = N by substituting the closed expression ()24II27|) for 
the transfer matrix elements into the definition (fTK|) and (|17)1. Exact results for \t\ 2 and |r| 2 
for the perfect system with absorption or gain, for E = 0, have been discussed by Datta^ 
and more extensive numerical results which include the additional effect of a weak disorder 
on the averaged logarithmic transmittance have been presented by Jiang and Soukoulis^. 
Special attention has been paid i n 12 ' 13 to the domain of intermediate lengths (in particular 
the critical length L c ) where the transmittance of an amplifying system changes from an 
initial growth at short lengths to an exponential decay at long lengths. 

In the following we discuss detailed results for transmittance and reflectance in the frame- 
work of the general analytic treatment for weak disorder in Sect. II. We shall consider suc- 
cessively the short- and the long lengths domains defined below. Our consistent treatment 
of the effect of a weak disorder in the framework of an exact analysis of absorption or gain 
at zeroth order leads to the identification of the important effects induced by absorption or 
amplification in the statistics of wave transport. 
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A. Short lengths 



For a fixed magnitude of the absorption/amplification parameter b the short lengths 
domain is defined by 

N\b\ « 1 , (45) 

or, equivalently L << l where l = 1/b for b > is the absorption length (in units of a) 
and Zo = —1/b for b < is the amplification length. We wish to obtain the logarithmic 
transmittance in the limit (|45|). which determines the short-length localization length. For 
this purpose we use the following approximations of (|4(JH42j) valid to lowest order, for small 
|6| and small m\b\, 

A m = i m (l-mb) , (46) 
D m = (-i) m (l + mb) , (47) 

B m = C m = h -{i m + {-i) m ) . (48) 

Note that these expressions would not be sufficient for discussing the reflection coefficients 
whose explicit forms differ for even and odd N and require inclusion of higher orders in 
ma = m for B m and C m . For brievety's sake we omit discussing the short lengths reflection 
coefficients in more detail. We substitute ()46ll48jl in the expression f!32[) for E = 0, which 
we then insert in fJHfi)) . After averaging over the disorder, using ()18j) . and performing the 
remaining geometric sums over sites, we obtain the following final results: 

(ln|t| 2 ) = -2bN- ^(l-Ab)N + 0(b 2 N 2 ) , (49) 
The short-length localization length obtained from ()35|) and ()49|) namely, 

1 , er 

yields 



- = 6 + ^(1 -46) , (50) 



7 = r + T-7T ' (51) 
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for absorption, and 



1 



1 1 4 



(52) 







for amplification, where 




(53) 



o 



Note that ()53|) is the exact perturbation expression (for E = 0) of the localization length 
for weak disorder, for 6 = 0. Indeed it coincides with the well-known exact result, £o = 
96W~ 2 sin 2 k, obtained by Thouless 2 ^, if the variance W 2 /12 of the rectangular distribution 
of width W of site energies h)M> is replaced by the gaussian mean square e^. The first 
two terms in (j51H52J) agree with the form of the short-lengths localization lengths derived 
previously from invariant imbedding**^. 

Finally, it is useful to clarify the general meaning of short- and long-lengths localization 
lengths in the framework of our weak disorder analysis. We recall that the perturbation 
treatment of disorder in Sect. II is valid for 



This implies, in particular, that the localization length ()53|) is valid for values el —>■ such 
that the limit (|5*4*|) embraces asymptotically large N for which the localization length is de- 
fined in ()35jl . Similarly, the short-lengths localization length (|5()j) in the presence of absorp- 
tion or gain is a true localization length only if it corresponds to the limit of asymptotically 
large N in (jHSJ). Thus if |fe| < el this limit is obtained for \b\ — > (since iV|6| < Nel « 1) 
while if |6| > el it is obtained by letting el — > 0. Now, for |6| < el (j^U) automatically implies 
(145)1 . in which case the short-lengths expressions ()51II52)) give the true localization lengths. 
On the other hand, for |6| > £q t wo possibilities exist for the localization lengths: 

• if for asymptotic lengths obeying (JSUK^o — * 0) one a l so nas ^1^1 <<: 1) then the 
localization lengths are clearly given by the "short-lengths" expressions (j51M52)l . This 
happens for values |6| > el sufficiently close to £q- 

• if for asymptotic lengths ()54|) the long lengths condition 



Nel « 1 



(54) 



N\b\ » 1 
12 



(55) 



is fulfilled, then the localization lengths are given by the equations ()64.al64.b|) and their 
limiting forms ()66.al 166. b|) and (j68.al68.bj) in subsection B below. This situation exists for 
|6|-values sufficiently larger than ejj. 



B. Long lengths 



The transfer matrix elements (j4(W42|) of a perfect system depend on the imaginary expo- 
nentials e ±iqN , which for |6| << 1 are given by 



D ±iqN _ / ±i xJV e =FJV(6+V+-) 



(56) 



where e iqN is exponentially growing for b < (amplification) and e~ tqN is growing for b > 
(absorption) in the long lengths regime, N >> (f53j). We first discuss the detailed 

form of the logarithmic transmission coefficient, ln|t| 2 = In t + c.c, and of the reflection 
coefficient |r ±T | 2 given by (fTTjl . (|TH|) and (j3()H34j) in terms of the transfer matrix elements 
An, Bn,Cn, Dn in (j4()H42|) . Retaining only the leading exponential terms at long lengths 
for absorption and amplification, respectively, we obtain successively 



for both signs of b, and 



In \t\ 



-2 



b 3 



N + 0(N 



(57) 



— ,b> 



(58) 



b 2 



b<0 



(59) 



The main feature of these results is that In \ t\ 2 is decreasing at large L for amplification as 
well as for absorption in agreement with previous studies^iiii2iii. 

Next we consider the effect of weak disorder at E = in the mean logarithmic transport 
coefficients (j36H38|) involving zeroth- and first order transfer matrix elements defined in ()401 
WB) and (|31H34|) . Using ()18|) , the averages of the various quadratic forms in first order transfer 
matrix elements in (|36H38|) reduce to simple sums over lattice sites m of products of two 
terms of the form M m _i — N m _i corresponding to the site m — 1 multiplied by a product 
of two terms of the form P/v_ m — Rn-tu corresponding to the site N — m (with M, N, P, R 
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representing elements, distinct or not, of the set of transfer matrix elements, A, B,C, D, 
of the pure system). Using (j40H42|) . we approximate the mth term in a given sum by the 
contribution which is independent of m, which yields the leading effect proportional to N 
for any of the sums involved (the contributions ignored in this approximation are readily 
shown to be of relative order -^). For the averages of products of first order transfer matrix 
elements entering in (|36H38|) we thus obtain the following results valid at E = 0, for any sign 
of b: 



((Q22) 2 ) 



Ne 



[{u + + vfe 



-iq(N-l) 



(60) 



tt_ — v e 



-iq(N-l)] 



(61) 



((QliY) = ((QnY) 
From dS2P and (JSZIEHD, it follows that 



(62) 



(ln|r+-| 2 ) = (ln|r-+| 2 ) = (ln|rf) 



(63) 



for both signs of b. 

Next we insert (j6()H62|) . together with (j41H42|l in (jrS5H38j) and simplify the resulting ex- 
pressions by retaining in each one of them only the leading exponential terms for iV|6| >> 1, 
successively for b > and b < 0. In this way we obtain the following exact expressions valid 
for arbitrary \b\ larger than and such that iV|6| >> 1: 



1 



2N 



In \e~ iqN \ 2 — £ o ( M + + ^j_ e 2 »g 



u 



+ 



(64.a) 



(In |r| 2 ) = In 



u + 



2 elN 



[U- — V 



e 2iq , b>0 ; (65.a) 



1 1 



ln\e iqN 



2 4 ( u - + v) 1 



-2iq 



U_ 



(64.b) 
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(In |r| 2 ) = In 



4 N / ^2 



(u + — v) 2 (u_ + v)' 



e~ 2iq , b < 



(65.b) 



For weak absorption/amplification |6| << 1 (with, however, iV|6| >> 1) we expand 
fl64.al65.aj) and ()64.bl65.b|) in powers of b, using (|43H44|) and (JSTIj) . To order b 2 in the effects 
of the disorder we obtain finally 



- = b-\ h — (1 

£ 3 8 V 



(66.a) 



, h 2 \ b 2 f 2 b 2 



(m|^) = m( T -- + 



(67.a) 



f 
1 



b 3 



+ 



(66.b) 



(ln|r| 2 ) = ln(^+| ; b<0 



(67.b) 



On the other hand, for strong absorption/amplification, |6| >> 1, we obtain successively 
from (j64.al65.aj) and ()64.bl65.b|) . to the orders indicated, 



1 . 6 3 el 



(68.a) 



(ln|r| 2 ) 



2 / J_ 

'6 V + 66 2 



26 2 



1 



6 > 



(69.a) 



1 



- 6 + 



+ 



86 2 



(68.b) 



(ln|r| 2 ) 



2 / 1 
6 V 1 ~ 3V 2 



+ 



2b 2 



1 + 



b<0 



(69.b) 



using expansions of (|43M4|) in powers if m. 

Our detailed results (j66.al66.bj) and ()68.al68.b|) for localization lengths and ()67.al67.b|) 

and ()69.al69.b|) for logarithmic reflection coefficients display remarkable new features related 
to the effects induced by absorption/amplification in statistical averages over the disorder. 
We recall that our results are valid at asymptotic lengths for values |6| > but not too close 
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to Sq (see the discussions in III. A above). In the absence of the induced statistical effects 
the results for inverse localization lengths coincide with the previously known results (1), 
(3)&*- On the other hand, in the absence of disorder, the results for the localization length 
(transmission coefficient) and for the reflection coefficient coincide with the exact results 
obtained by Datta^. In particular, in the absence of disorder the reflection coefficient is 
asymptotically constant for N — > oo 12 i 13 . 

Concerning now the statistical effects induced by absorption/amplification in the inverse 
localization lengths, our results ()66.al66.b|) and ()68.al68.pp lead to the following conclusions: 

1. the effects are identical for absorption and for amplification, for weak- as well as for 
strong absorption/amplification. 

2. the statistical effect induced by absorption/amplification increases the localiza- 
tion length for weak absorption/amplification while reducing it for strong absorp- 
tion/amplification. 

3. localization by disorder is destroyed in the presence of sufficiently strong absorp- 
tion/amplification. 

On the other hand, our results in (|67.al67.b|) and (|69.al69.b|) for the statistical effects 
induced by absorption/amplification in the reflection coefficient reveal that: 

1. a weak absorption induces a weak asymptotic statistical growth of (\r\ 2 ) while a cor- 
responding weak amplification leads to no statistical effect. 

2. for large absorption/amplification parameters |6|, absorption and amplification induce 
identical weak statistical growth terms (to leading order in \b\~ 2 ) in (|r| 2 ). 

IV. CONCLUDING REMARKS 

The main results of this paper are summarized in the analytical expressions (j51H52j) 
()66.al l66.bM69.al69.bj) for inverse localization lengths and logarithmic reflection coefficients 
in short- and long random tight-binding systems with absorption or gain. These results are 
discussed in detail in the main text. Our analysis in Section II is valid for Neq << 1 which 
characterizes the weak localization regime identified more generally by the limit L << £ - 
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It would be interesting in future work to study the effect of absorption or amplification in 
the strong localization (or localized-) regime L » £o, for weak disorder. This would allow 
to study the additional effects associated with anomalies in £ existing at special energies, 
in particular at the band center—. The study of transmission and reflection in the localized 
regime requires a more involved treatment of the disorder, respecting, in particular, the 
asymptotic unitarity limit of the reflection coefficient in the absence of absorption or gain. 
A simple analytic treatment of statistical properties of the transmittance in the localized 
regime in the absence of absorption has been discussed recently in2i. 

We close with brief remarks on the respective roles of different symmetries of the transfer 
matrix (or of the lack thereof) for the disordered tight-binding system with absorption or 
gain studied above. In Sect. II a central role is played by the transfer matrix Q n of an 
elementary disordered segment enclosing just one site n. Q n obeys the property 

detQ n = l , (70) 

which leads to the relations 

U++|2 = u—|2 I +-I2 I -+I2 (n) 

for the reflection- and transmission coefficients for waves incident from the left and from the 
right, respectively. As shown in Sect. II, (fTUj) implies that det Q = 1 which in turn leads to 
the identity of the transmission coefficients |t ++ | 2 and \t | 2 (Eq. 1)17])) for a system of N 
sites. Now (J71|) may be viewed simply as reflecting the symmetry of the piecewise defined 
solutions of the Schrodinger equation for plane waves incident from the right and from 
the left, respectively, in a single-site random segment. This finally shows that the transfer 
matrix Q embodies the basic left-right symmetry of Eq.(j2J), via Eq. (jl6|) . which leads to 
the properties (\t ++ \ 2 } = (\t | 2 ) = (|t| 2 ) (|r + ~| 2 ) = (|r~ + | 2 ) = (|r| 2 ) for the observable 
transmission and reflection coefficients. 

In the absence of absorption or gain the disordered system (J2J) possesses a further well- 
known symmetry, namely time-reversal symmetry. This symmetry implies that the 2x2 
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transfer matrix X satisfies the condition 

- ~ (o A 

A" = aXa, <r = . (72) 

V 10 / 

This symmetry is broken when absorption or gain is present as follows e.g. from the transfer 
matrix (|23|) for the pure system. Indeed, for E — 0, (J72*j) would require that 

B* N = C N and A* N = D N , (73) 

which is not the case for the elements (j40H42J) . Now it is readily seen that the lack of time- 
reversal symmetry as shown by the violation of ()73)1 is related to a physical fact, namely the 
absence of current conservation, which means that \r\ 2 + \t\ 2 ^ 1. Indeed, from (|17jl and 
(|2~n|) we have in the present case 



>f + w 2 = 1 4¥i L * 1 • («) 



as seen from (|2*5jl . since A^r 7^ and ^ B* N . Note also a further related consequence 
of the lack of time-reversal symmetry of the perfectly absorbing or amplifying systems: this 
is the violation of the duality relation for the scattering-matrix derived by Paaschens et al. 9 . 
Violation of the duality relation of^ for the 5*-matrix is easily demonstrated by substituting 
the transfer matrix elements (|41H43|) in Eqs. f)14H15j) for the reflection and transmission 
amplitude coefficients. 
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